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II. NOMENCLATURE 




D 

L 

M 

00 

P 

q. 

r 

o 

*1 

R 

c 

U 

U 

CO 

V 

X 

Y 
a 

U) 

v 



Magnus force coefficient based upon frees tream velocity and 

side area 

2C y /5p 

Sideforce correlation constants for small and large at 

respectively 

Cylinder diameter 

Cylinder length 

Freestream Mach number 

Dimensionless spin rate (v/u) 

Freestream dynamic pressure 
Cylinder radius 

Reynolds number based on cylinder length 

Crossflow Reynolds number based on crossflow velocity and 
diameter 

Freestream velocity of spinning missile 
Freestream velocity of rotating cylinder at 90° ot 

Cylinder surface speed 

Longitudinal distance from nose 

Sideforce 

Angle of attack 

Spin angular velocity 

Kinematic viscosity 
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III. INTRODUCTION 



The production of sideforce (Magnus effect) on spinning bodies at 
angle of attack has been of interest to aerodynamicists since the l800's 
when G. Magnus published the first experimental results of the drift of 
musket shot due to spin. Since that time, it has been recognized that 
accurate prediction of Magnus forces was necessary to calculate the tra- 
jectory and stability of spinning shell and missiles. In spite of all the 
time and effort expended since the l800's, theoretical predictions of 
Magnus force production are extremely crude. Although the basic mechanism 
for Magnus force production is recognized, the mathematical model required 
for solutions considering bodies with even reasonably simple shapes is too 
complex and time consuming even for today's sophisticated fluid mechanics 
calculation methods. The engineer is forced, then, to use the very simple 
phenomenological models developed to date in conjunction with experimental 
correlations to make estimates of the Magnus effect on spinning missiles. 
This report gives the details of such a study. 
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IV. MAGNUS FORCE PREDICTIONS 



It is generally agreed that Magnus force production is due entirely to 
viscous effects. An asymmetrical boundary layer is produced on the missile 
surface by the combined effects of spin and angle of attack. At very low 
angles of attack the boundary layer is attached over most of the missile 
length so that sideforce predictions may be based upon displacement thick- 
ness and shearing stress calculations. At larger angles of attack the 
asymmetrical boundary layer separates causing a shedding of vorticity into 
the outer inviscid flow. The shed vorticity induces an equal but opposite 
circulation in the outer flow producing an additional circulation sideforce 
contribution. The modeling of such a three-dimensional boundary layer 
separation problem is beyond the capabilities of present day theoretical 
techniques even for the relatively simple case of an incompressible, 
laminar flow over a rotating body of revolution with a cylindrical after- 
body. Theoretical predictions of the effect of compressibility and a tur- 
bulent boundary layer on Magnus force must necessarily be even more crude. 

At low angles of attack and small nondimens ional spin rates the most 
realistic incompressible, laminar flow theory to date is that of Kelly."*" 
Kelly considered the outer flow of an open ended rotating circular cylin- 
der at an angle of attack. A perturbation solution of the boundary layer 
equations was used to calculate displacement thickness and shearing stress 
contributions to the Magnus force. The effects of a forebody and shed 
boundary layer vorticity were not considered. At very low angles of 
attack and for a high fineness ratio body this linear theory yeilds a rea- 
sonable estimate of the Magnus force. 

For larger angles of attack (but still small) Power and Iverson 
made an estimate of the Magnus force contribution due to boundary layer 
shed vorticity. This theory will be outlined in more detail because the 
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sideforce experimental data correlation used in this report originated 
from this work. 

Following Kelly's lead, the incompressible, laminar flow about an 
open-ended circular cylinder at angle of attack shown in figure V-l was 
considered. It was assumed that the effect of shed vorticity could be 
modeled by using a crossflow analogy. The circulation induced into the 
outer flow was approximated by assuming that the crossflow along the longi- 
tudinal cylinder axis was similar to the unsteady solution for the flow 
over an impulsively rotated cylinder placed perpendicular to the free- 
stream as shown in Figure V-2. Values of the time dependent lift coeffi- 
cient for the impulsively rotated cylinder where estimated from the results 

4 

of Thoman and Szewczyk. It was found that this data correlates very well 
by the relationship: 




(IV. 1) 



where U ffl is the frees tream velocity, D is the cylinder diameter, 
the cylinder surface speed. The crossflow analogy assumes that: 

U = U sin a 



U cos a 

_L = _L_ 

U sin a 



and V is 



(IV. 2) 



Using the additional assumption that f(jj~) can be approximated by a linear 

00 

. function the local Magnus sideforce coefficient based on diameter and 
crossflow velocity is given by: 



C y (x) 



= C 





x JL 


/ - \ 
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(IV. 3) 
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where C is a constant, 
o 

For small angles of attack equation (V.3) leads to an equivalent local 
circulation strength given by 



k(x) 




(IV. 4) 



where is a constant, R c is the crossflow Reynolds number based on cross- 
flow velocity and the cylinder diameter, and T o is the surface circulation 
strength. The total Magnus sideforce was calculated by integrating the 
local circulation contributions over the length and adding the contribution 
due to displacement thickness and shearing stress asymmetry. The result 
for small a and nondimens ional spin rate p = V/U is: 

C y = k^ + k 2 c 7 /4 (IV. 5) 



where : 
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and the ffegnus coefficient is based on the side area LD. 

Equation (V.5) predicts that for a given missile configuration, spin 
rate, speed, and altitude, the sideforce coefficient is nonlinear with <y. 
The constant can be evaluated by considering experimental data for 
cases which are consistent with the assumptions required for crossflow 
analogy application. It would also be expected that might be a func- 
tion of Mach number because of compressibility effects. Because of the 
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low crossflow velocities at low angles of attack, crossflow Reynolds num- 
bers are low enough for the laminar assumption to be valid even at high 
forward speeds. 

Equation's (V.5) nonlinearity suggests that the dominant Magnus force 
production mechanism is shed vorticity at higher angle of attack. If 

5 

this assumption is made, a much simpler correlation suggested by Iverson 
can be developed. Returning to the result of equation V.3> an integration 
of the circulation produced local sideforce results in a Magnus force 
given by: 

C Y = C q p sin a 

where 

tan ot /y. ]t • ii 

u = — — — /R 4 sm 4 a 
D ' d 

By expanding the exponential we find that 



1 / -uL , 

1 + 5I (e - V 
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For reasonable Reynolds numbers and high fineness ratios the Magnus force 
may be approximated by: 
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tan a sin 
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(iv. 6) 



Equation V.6 predicts that all Magnus data should collapse along a line of 
slope C q if SC^/dp = is plotted against the correlation parameter 



1 , .3/4 /l\ 
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FIGURE IV. 2: SPINNING CYLINDER CROSSFLOW ANALOGY 



V. EXPERIMENTAL MAGNUS MEASUREMENTS 

Theoretical predictions of Magnus side force based upon vorticity 
shed into the outer flow show nonlinear variation with angle of attack. 
Experimental Magnus measurements were used to compare actual side force 
production to that predicted by theory and to establish the values of 
the empirical factors in equations V.5 and V.6. For this purpose the 
model configuration selected for study had a tangent-ogive nose followed 
by a constant 2.95 inch diameter cylindrical afterbody. The afterbody 
length was varied to provide three missile length-to-diameter ratios, 
l/d, of 9*9> 7.8 and 6.0. 

The model was machined from aluminum and carefully balanced to re- 
duce vibrational loads when spinning. Figure VI -1 shows schematically 
the model, its support system, the balance, and the air motor drive 
system. The l/8 horsepower Standard Pneumatic rotary- vane, bi-directional 
compressed air motor allowed missle rotational rates as high as 3000 RPM 
to be tested. An electromagnetic tachometer pickup was used in conjunc- 
tion with a Monsanto 101B counter to determine model spin rate. 

A simple beam balance was used to measure the Magnus side force and 
center of pressure location. Two four-element strain gage bridge circuits 
were mounted on the beam flexure such that only sideforce and Magnus 
moment were measured. Air motor pressure lines were in the vertical 
plane of the balance to reduce pressurization effects on balance output. 

A subsequent balance calibration showed that for pressures less than 80 
psia the interaction of line pressure and sideforce measurement was 
negligible. Figure VT. 2 shows the wiring diagram for the balance system. 
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The bridge circuit outputs were amplified and displayed on a HP 7100 B 
two channel chart recorder. The sting balance system was calibrated by 
hanging known weights at known distances from the strain gage locations 
and plotting the applied moment verses chart recorder output. The data 
reduction program in the Appendix was used to calculate the aerodynamic 
coefficients. Figure VI-3 shows the model mounted in the wind tunnel. 

Data was taken at a constant nominal Reynolds number based on missile 
length of 1 x 10 and at angles of attack up to 15 degrees. Sideforce 
production at seven nondimensional spin rates between +0.4 was measured 
at each angle of attack. Air motor supply pressure was adjusted until 
the desired spin rate was established. Figure VI-4 through VI-6 show 
the resulting Magnus force coefficients and the location of the Magnus 
force center of pressure with respect to the missile base. These figures 
show a nonlinear variation of sideforce coefficient with angle of attack 
at constant nondimensional spin rate. The Magnus center of pressure is 
seen to be approximately one third of the missile length from the base for 
all the cases tested. Center of pressure location was insensitive to 
angle of attack and spin rate. The L/D =9-9 data, however, show movement 
of the center of pressure with both angle of attack and spin rate for 
angles of attack less than ten degrees. Sideforce coefficient at constant 
angle of attack and spin rate is dependent upon missile L/D. Figure VI-7 
shows this effect of L/D on the Magnus sideforce coefficient for a non- 
dimensional spin rate of 0.2. 

Theory predicts a linear behavior of sideforce coefficient with 
nondimensional spin rate p. Figure VI-8 presents the experimental 
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data plotted in such a way to test this result. The experimental side- 

force coefficients were divided by the nondimens ional spin rate to 

calculate C^_ = dCy/^p. If sideforce coefficient is linear in p all the 
P 

data should collapse to a single line when C y _ is plotted against angle 

P 

of attack. Figure VI-8 shows that for a constant missile L/D the data 
show an excellent agreement with this theoretical result. Also shown 
on Figure VT-8 is the linear theory prediction of Kelly'*’. The data 
show excellent agreement with the linear theory at angle of attack less 
than five degrees. At angles greater than five degrees the data show 
nonlinear behavior with angle of attack. 

An attempt to correlate the data with equation V-5 was unsuccessful. 
The data reduction program calculated the value of the correlation constant 
for each experimental point. Values of ranging from zero to three 
resulted. The value of increased as missle L/D and angle of attack 
increased. No single value of C fits the experimental data. Figure 
VI-9 shows a comparison of equation V-5 with the value of taken as 
unity. Equation V-5 consistently overestimates the Magnus coefficient at 
angles of attack less than approximately ten degrees and underestimates 
the Magnus force at higher angles. 

Figure VI-10 shows the experimental data plotted against the 
correlation parameter of equation V-6. Excellent agreement with this 
correlation was found. The data of the higher L.D ratio missle shapes 
collapsed about a straight line of slope 6.6. The L/D = 6.0 data also 
collapsed as predicted but at a value of the correlation constant equal 
to six. 



13 




I— I 

<u 



£ 



-P 

w 

<u 

H 

a) 



cCj 

-P 

<h 

o 



•rH 

cd 

-p 

a) 

Q 



> 

a) 



•H 



14 




15 



Figure V.2: Wiring Schematic 
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Figure V.3: L/d = 9-9 Test Model 




SIDE FORCE COEFFICIENT - C v CENTER OF PRESSURE LOCATION 





Figure V.4: Magnus Force and Center of Pressure for a 9*9 L/D 

Spinning Missile 
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Figure V.5: Magnus Force and Center of Pressure for a 7.8 L/D 

Spinning Missile 
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SIDE FORCE COEFFICIENT - C v CENTER OF PRESSURE LOCATION 



0.6 



0.4 - 




0.2 - 




0 



- L - l L 

5 10 15 



ANGLE OF ATTACK - a (degrees) 




ANGLE OF ATTACK - a (degrees) 



Figure V.6: Magnus Force and Center of Pressure for a 6.0 L/D 

Spinning Missile 
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MAGNUS COEFFICIENT 




Figure V.8: Magnus Force Coefficient C^_ 

P 



21 




MAGNUS COEFFICIENT 




ANGLE OF ATTACK - 0/ (degrees) 



Figure V.9: Magnus Force Comparison with Equation V.5. 
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0.2 L R = 0.967 X 10 
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VI. CONCLUSIONS 



A comparison of the experimental data in the previous section 
with theory leads the author to the following conclusions: 

1. For the configurations tested the Magnus coefficient 
is linear in spin rate for constant R , <y and L/D for nondimensional 
spin rates less than +0.4. 

2. Kelly T s linear theory can predict accurately the 

variation of C^ with l/D and only at very low angles of attack 
P 

(in the case of the configuration tested the angle of attack must 
be less than five degrees.) 

3. The crossflow analogy prediction that Magnus coefficient 
is nonlinear with angle of attack is accurate. The data, however, 
only match the theory if the linear term in angle of attack is 
neglected. It might, therefore, be concluded that although the 
Magnus effects due to the attached boundary layer growth (displacement 
and shearing stress contributions) are important at very small angles 
of attack, the primary mechanism of Magnus force production is through 
shed vorticity even at small angles of attack. 

4. The crossflow analogy leads to an accurate Magnus force 
correlation parameter if L/D is sufficiently high (in the case of the 
present data L/D must be greater than approximately seven). 

5. Iverson^ has shown that the Magnus force coefficient for 
supersonic freestream Mach numbers also show good correlation until 
the crossflow is near sonic. The correlation constant C q is expected 
to be an increasing function of the crossflow Mach number until this 
velocity is reached. At this point the- parameter C q will remain relatively 
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constant because of the inability of shed vorticity to reach the outer 



supersonic flow. The experimental value of C q for incompressible 
crossflow Mach numbers has been determined by this study 
to be 6.6 if the coefficient is based upon side area. The experimental 
data therefore correlate well with the prediction: 

C Y _ = 3.3 (L/D)5/yA 



for angles of attack less than fifteen degrees. 
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COMPUTER PROGRAM FOR DATA REDUCTION 



READ (5,100) B,TCF,RK1,RK2,N 
100 FORMAT (5E10.0) 

READ (5,100) X,D,PA,TF ,R1L 

C B IS DISTANCE BETWEEN GAGES (IN) ,TCF IS THE TUNNEL CALIBRATION 

C FACTOR (CM H20/PSF),K1 & K2 ARE BALANCE CALIBRATION SLOPES 

C (IN LB/ CHART LINES ),N IS NUMBER OF DATA POINTS 

C R1L IS DISTANCE FROM FRONT GAGE TO BASE (IN) 

C X IS MISSLE LENGTH (FT),D IS DIAMETER (FT), PA IS ATMOSPHERIC 

C PRESSURE (PSF) ,TF IS TEST SECTION TEMP.(F) 

C Cl IS SMALL ALPHA CORRELATION CONSTANT AND C2 IS FOR LARGE ALPHA 

DO 1 1=1, N 

READ (5,100) R1,R2 ,RPM ,A,DP 

C R1 IS FRONT GAGE READING (CHART LINES ) ,R2 IS BACK GAGE READING 

C RPM IS ROTATIONAL SFEED (REV/MIN.),A IS ANGLE OF ATTACK, DP IS 

C DELTA TEST SECTION PRESSURE (CM H20) 

Y= ( RK2*R2 -RK1*R 1 ) /B 
XCP=RKl*Rl/Y 
XCPR=(XCP+R1L) / (X*12. ) 

C Y IS SIDE FORCE (LB), XCP IS CENTER OF PRESSURE LOCATION (IN) 

C XCPR IS NONDIM CP LOCATION FROM BASE 

TR= TF +460. 

RHO= PA* 70 . 73/ ( 1716 . *TR ) 

ESB= 0.01295* (0.868 +.4*D/x)*D*D*X 
Q=DP* ( 1 . + 2 . *ESB ) /TCF 
V= (2.*Q/RH0)**0.5 

C CY IS SIDE FORCE COEFF BASED ON SIDE AREA 

CY=Y/Q*D*X) 

PHAT=3 . l4l6*D*RPM/ ( 60 . *V ) 

xmu =( 2 . 26967 *tr** 1 . 5 )/ ( (TR+ 198 . 72 )* 10 .** 8 . ) 

RL=RHO*V*X/XMU 
XLX= X/D 
RAD= A/57.3 
R0= RL*SIN(RAD) /XLX 

XK3= XLX*(60. 97-357. 37*XLX/(RL**. 5) -227. 01*XLX*XLX/RL)/ 
l(RL**.5)*3.l4l6/4. 
sa=sin(rad) 
ta=tan(rad) 

C CYCP IS THE C ORREA LAT ION PARAMETER FOR LARGE ALPHA 

CYCP=XLX**1 . 25*SA** . 75*TA / ( 2 . *RL** . 2 5 ) 

IF (PHAT) 3,2,3 

2 CYP=0 . 0 
C1=0.0 
C2=0 . 0 
GO TO 4 

3 CYP=CY/PHAT 

IF (RAD) 6,5,6 • 

6 IF (SA) 7,5,7 

7 IF (TA) 8,5,8 
5 C1=0.0 

C2=0 . 0 
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GO TO 4 

8 Cl=4 . *RL** . 25* ( CY/FHAT -XK3*RAD )/ ( 3 . l4l6*XLX**l . 25*RAD**1 . 75 ) 
C2=2 . *CY8RL** . 25/ (PHAT*SA** . 75*TA*XEX**1 . 25 ) 

4 WRITE (6,110) 

110 FORMAT (IH,// 19X, 'ALPHA' ,16X,’PHAT' ,19X,'CY' ,17X,'XCPR' ,18X, 
1 ' CYP ' , 17X , ’ CYCP ' ) 

WRITE (6,111) A,PHAT,CY,XCPR,CYP,CYCP 

111 FORMAT ( 9X,6E20.7) 

WRITE (6,112) 

112 FORMAT (IH,/ 39X, ' L/D ’ ,l 8 X, 'RL' ,l 8 X, 'RC ’ ,17X, ' Cl ' ,l 8 X, ' C2 ’ ) 
WRITE (6,113) XLX ,RL ,RC , Cl , C2 

113 FORMAT (29X,5E20.7) 

1 CONTINUE 

STOP 

END 
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